Limit cycle theory of temporal current self-oscillations in sequential tunneling of 

superlattices 
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A unified theory of the temporal current self-oscillations 
is presented. We establish these oscillations as the manifes- 
tations of limit cycles, around unstable steady-state solutions 
caused by the negative differential conductance. This theory 
implies that both the generation and the motion of an electric- 
field domain boundary are universal in the sense that they do 
not depend on the initial conditions. Under an extra weak ac 
bias with a frequency uj a c, the frequency must be either ui ac or 
an integer fractional of uj ac if the tunneling current oscillates 
periodically in time, indicating the periodic doubling for this 
non- linear dynamical system. 
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Following the discovery of temporal current self- 
oscillations (TCSOs) in sequential tunneling of superlat- 
tices (SLs) under a dc bias [1-3] , a large number of exper- 
imental and theoretical studies have focused on their ori- 
gin and how these oscillations develop from steady-state 
solutions (SSSs). Experimentally, TCSOs have been ob- 
served in both doped and undoped SLs [1-3] . The oscil- 
lations can be induced by continuous illumination with 
laser light [3] or by a change in doping [1]. Recently, it 
has been shown that TCSOs can also be induced by ap- 
plying an external magnetic field parallel to SL layers or 
varying the sample temperature [4] . 

Our current theoretical understanding of TCSOs is 
mainly from numerical studies [5] . Early works of Bonilla 
and his co-workers [5] established a correct model for TC- 
SOs. Many numerical results and some analyses were 
also done. They simulated and reproduced many exper- 
imental results, including finding proper model parame- 
ters to simulate different experimental situations. Great 
progress in understanding of TCSOs had been made be- 
cause of them. However, a simple physical picture did 
not appear in those early works. The understanding at 
the computational level is the first step, and deep in- 
sights can be only obtained when the general concepts 
and principles are found. There are also microscopic 
Green's function calculations [6]. While the microscopic 
approach would be accurate if all the microscopic pa- 
rameters and mechanisms were known, it remains a chal- 
lenge to deduce the rules of macroscopic behavior from 
the microscopic details. Recently, a clear route to TC- 
SOs developed from SSSs in sequential tunneling of SLs 
was proposed [7]: Due to the negative differential con- 



ductance (NDC) [8], a SSS is not stable. A limit cycle is 
generated around the unstable SSS because of the local 
repulsion and global attraction in the phase space. The 
system moves along the limit cycle, leading to a TCSO. 

Unfortunately, more careful studies [9,10] show that 
the simple model used in reference 7 does not have a 
TCSO solution even though it gives the on-set of insta- 
bility conditions of a SSS. Thus, the evidence for limit 
cycles in the TCSO regime is still lacking. In this letter 
we search for the fundamental concepts and general prin- 
ciples for TCSOs, based on the correct model of Bonilla 
and his co-workers [5]. We apply the concept of limit 
cycle, and explain TCSOs as the manifestations of limit 
cycles. Although numerical results are not our emphases 
in this paper, we show that they can be easily under- 
stood by using this concept. For example, the generation 
and motion of an EFD boundary do not depend on the 
initial conditions. An EFD-boundary does not necessar- 
ily start from the emitter, and end up in the collector. 
Furthermore, limit cycle concept gives us a power of pre- 
diction. We predict that the generation and motion of 
EFD boundaries are universal, and the frequency of a 
periodic motion under an ac bias must be either the ac 
bias frequency or its integer fractional. The limit cycle 
is a well-known concept in non linear physics. Thus, we 
find that TCSOs can be understood under the general 
concepts of non linear physics. 
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FIG. 1. Schematic illustration of an SL system, fn is the 
local chemical potential of the i th well, /Ltz, and (m are the 
chemical potential of left-hand side and right-hand side elec- 
trodes, respectively. Vi is the bias on the ith barrier, and 
fiL — jJ.R = U is the external bias. 

We consider a system consisting of N quantum wells 
as shown schematically in Fig. 1. An external bias U is 
applied between the two end wells. Current flows per- 
pendicular to the SL layers. In the sequential tunneling, 
charge carriers are in local equilibrium within each well, 
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so that a chemical potential can be defined locally. The 
chemical potential difference between two adjacent wells 
is called bias V on the barrier between the two wells. A 
current Ii passes through the i th barrier under a given 
bias Vi. This current may depend on other parameters, 
such as doping Nd. One of the results in reference 7 is 
that a SSS must be unstable if there are two or more 
barriers being in the NDC regime. It is worth pointing 
out that, although this instability result is obtained from 
an analysis of a simplified sequential tunneling model, it 
is generally true. Without losing generality, we assume 
that only barriers 1 and 2, which create well 2 shown in 
Fig. 1, are under a NDC regime for a SSS. Assume the 
chemical potential in well 2 increases a little bit due to a 
fluctuation. Then bias V\ on barrier 1 decreases while Vi 
increases. Because both of the barriers arc under NDC 
regime, charge carriers flow more into well 2 through bar- 
rier 1 while less carriers flow out of it, leading to a further 
increase of the chemical potential in well 2. This drives 
the system away from the SSS, i.e. instability. 

Following reference 5, the dynamics of the system is 
governed by the discrete Poisson equations 

k(V i -V i - 1 )=n i -N D , i = l,2,...JV, (1) 

and the current continuity equations 
BV 

J = k-^+I h 2 = 0,1,2,. ..N (2) 

where k depends on the SL structure and its dielectric 
constant, rij is the electric charge in the i th well. In 
Eq. (1), a same doping in all wells is assumed. I, is, in 
general, a function of Vi and n,. It can be shown [9] that 
all SSSs are stable if Ii is a function of Vi only. On the 
other hand, a SSS may be unstable [5] if one chooses Ii = 
fiiv(Vi), where v is a phenomenological drifting velocity 
which is, for simplicity, assumed to be a function of Vi 
only. The constraint equation for Vi is 

J2Vi = U. (3) 

Previous studies [5] proved that this model is capable 
of describing TCSOs. To close the equations, a proper 
boundary condition is needed. It is proper to assume a 
constant no, tiq = SNd with S as a model parameter, 
if the carrier density in the emitter is much larger than 
those in wells, and its change due to a tiny tunneling cur- 
rent is negligible. However, it is mathematically equiva- 
lent to other boundary conditions used in literature [5]. 
Our goal is to show that the limit cycle is one of the most 
important features in this widely studied model. 

According to our theory, NDC is essential for TCSOs. 
Thus, a TCSO can only occur when there is a negative 
differential velocity in v(V) [5]. There are many ways of 
choosing it. One can assume v being a sum of a series 



of Lorentzian functions if NDC is due to the resonance 
tunneling between discrete quasi-localized states of wells. 
One may also choose a piecewise linear function in order 
to make an analytic investigation easy. We shall assume 
v(V) as the sum of two Lorentzian functions, v(V) — 
0.0081/[(V/.E - l) 2 + 0.01] + 0.36/ [(V/E - 2.35) 2 + 0.18]. 
This v has two peaks at V = E and V = 2.35E. A 
negative differential velocity exists between V = E and 
about V = I.3.E. Thus, E can be used as a natural 
unit of bias, and l/v(E) as that of the time (the lat- 
tice constant is set to be 1). For N = 40, U = 43.6E, 
Nd = 0.095kE, and S = 1.001, the set of equations has 
TCSO solutions. We solved the set of equations by using 
the Runge-Kutta method. Initially, the external bias is 
randomly distributed. Quickly, Vi reaches a stable state 
which does not depend on the intial conditions. Figure 2 
is the projection of the stable bias-trajectory in V5 — V38 
phase plane. V5 is in the low EFD while V38 is in the 
high EFD here. Clearly, one obtains a closed isolated 
curve indicating a limit cycle. The current oscillation 
period, the inverse of system intrinsic frequency uj , is 
the time that the system needs to move around the cycle 
once. The inset is the phase diagram of the TCSOs in 
U-Nd plane with all other parameters unchanged, where 
No is in units of kE. The shadowed area corresponds to 
the TCSO regime. Of course, the diagram depends on 
the values of 5, N, and function of v. 
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FIG. 2. The trajectory of the system in phase plane V5 — V38 
in a TCSO regime. The closed curve is the projection of the 
limit cycle in the plane. The bias is in units of E. The inset is 
the phase diagram in U-Nd plane with all other parameters 
unchanged, where Nd is in units of kE. The system inside 
the shadowed area will have a TCSO solution while it has a 
static current- voltage characteristic outside this area. 

Although it is known that TCSOs are accompanied by 
the motion of EFD boundaries, how EFD boundaries are 
generated and propagate inside SLs were debated [5] . Ac- 
cording to the limit cycle picture, the charge accumula- 
tion (depletion) is responsible to the creation of an EFD 
boundary. Charge carriers are accumulated (depleted) 
in a particular well because of an imbalance of carriers 
flowing in and flowing out. This imbalance is caused by 
NDC as we argued early. Thus an EFD boundary can 
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start at any well and oscillates inside a SL. Furthermore, 
the limit cycle picture means that EFD boundaries should 
not depend on the initial conditions, but are completely 
determined by the limit cycles around unstable SSSs. 

To demonstrate the correctness of our picture, we lo- 
cate numerically the position of the EFD boundary in 
the calculation that gives Fig. 2. Figure 3 is the evo- 
lution of the boundary. It reaches a stable oscillating 
state quickly. One can see that the EFD boundary os- 
cillating between wells 26 and 37 in the SL of total 40 
wells. The inset is the field distribution across the SL 
at points a, b, c, and d in Fig. 3. The bias V5 (V38) 
of the low (high) EFD moves up and down as the EFD 
boundary oscillates inside the SL. The stable oscillation 
state, a manifestation of a limit cycle, does not change 
when different initial conditions are used. In this sense 
both generation and motion of an EFD boundary are uni- 
versal. Although early numerical calculations [5] might 
have already implied that an EFD boundary can start in 
an interior well and oscillate inside a SL. It may not be 
easy for theories like that of reference 5 to explain this 
universal property. 
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FIG. 3. The time evolution of an EFD boundary. The 
parameters are the same as those for Fig. 2. The boundary is 
oscillating between well 26 and well 37. The stable oscillation 
does not depend on the initial condictions. The inset is the 
electric field distribution at points a, b, c, and d in Fig. 3. 

Except a few of attempts [11,12] which were mainly 
on the numerical aspects dealing with non-periodic time- 
dependent tunneling current, most theoretical studies 
have not considered TCSOs under the influence of an ac 
bias due to the lack of a clear physical picture such that 
one would not be able to analyze its effects. The limit 
cycle theory offers a way of analyzing the ac bias effect. 
Applying a small extra ac bias, it affects TCSOs through 
perturbing the system trajectories. In the case that the 
limit cycles still exist and are stable, or in other words, 
the tunneling current can oscillate periodically with fre- 
quency u, the ac bias on the system should return to its 
starting value after the system completes one-round mo- 
tion on the limit cycle. It means Lo ac /oj = n = integer. A 
weak ac bias cannot greatly change the evolution trajec- 
tory of the system in the phase space. The consequences 



are as follows: a) Current oscillation frequency uj cannot 
be much greater than the intrinsic frequency. Thus, at 
high ac bias frequency (w ac >> Wo), it is impossible for 
the bias to deform the limit cycle slightly such that the 
time for the system moving around the cycle once to be 
the same as the period of the ac bias, b) In a case that 
the system cannot deform itself to match io acl a trajec- 
tory may become a closed curve after several turns in 
the phase space. Therefore, uu /u = ra = integer. In- 
deed, figure 4 is the limit cycles in phase plane E^-E^ 
under an extra ac bias V ac cos(w ac i) with V ac = 0A4E 
and uo ac = 2uot) (dot line), 3^o (solid line) while the rest 
of parameters remain the same as those for Figs. 2- 
3. The right (left) inset is the current-time curve for 
u) ac = 2u; (3uj ) after the current oscillation becomes 
stable. The Fourier transformation shows the current os- 
cillation frequency u> being uq. This is the solution of 
dJ ac /uj = n = integer and ujq/uu = m = integer for the 
smallest possible n (=2, 3) and m (=1). 
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FIG. 4. The limit cycles in phase plane V5-V38 under an 
extra weak ac bias V ac cos(uj ac t) with Lo ac being 2ojo and 3o)o- 
The inset on the right (left) is the current-time curve for 
LOac = 2cjo (—3luo). V ac = 0A4E, and all other parameters 
are the same as those for Figs. 2 and 3. The time is in units 
of l/v(E). biases are in units of E. The current frequencies 
are io ac /2 = luq and L0 ac /3 = loo, respectively. 

For uj ac /uj being a non-integer rational number, u> 
should be different from uj according to the rules of 
uJacju = n = integer and ljq/oj = m = integer. For 
example, uj ac /u; = 1/q with q =integer, then the limit 
cycle can deform itself in such a way that it becomes 
a closed curve after q turns in the phase space, corre- 
sponding to n = 1 and m = q. In this case, the current 
oscillates with oj ac . The solid line in figure 5 is the nu- 
merical results of the limit cycle in phase plane V5-V38 
for u ac = uj /3 while the rest of parameters are kept 
the same as those for Figs. 2 to 4. Indeed, the limit 
cycle, which does not depend on the initial conditions, 
makes q = 3 turns in the phase plane as expected. The 
left inset is the corresponding tunneling current evolution 
curve. The Fourier transformation of the current evolu- 
tion indeed shows u = u ac . For u) ac /uu = 2.5, our rules 
predict the ui = ujo/2 = uj ac /5, corresponding to n = 5 
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and m = 2. This result is verified by the numerical cal- 
culation as displayed in figure 5 (dot line). As before, all 
other parameters are kept the same as those for Figs. 2 
to 4. The right inset is the corresponding tunneling cur- 
rent evolution curve. The limit cycle is a two-turn closed 
curve in phase plane V5-V38- The Fourier transformation 
of the current evolution demonstrates that the current 
oscillates with frequency wo/2 = u> ac /5. We would like 
to make following remarks, a) lo = co ac for co ac = uj /q 
cannot be true for all q because lo should approach loq 
in the limit u ac — > 0. b) We have considered only peri- 
odic responses of the tunneling current. It does not rule 
out other more complicated behaviors. In fact, there are 
reports [11-13] of chaotic current-time behavior under a 
combined dc and ac biases, c) Vary external parameters, 
the size and shape of a limit cycle should change in gen- 
eral. Thus lvq can shift. As a consequence, co ac /co = n 
may be quite robust against Lo ac . When oo ac is very close 
to loq, it may be possible for the limit cycle to deform 
itself in such a way that loq shifts to cj ac . In this case, 
the current shall oscillate with frequency uj ac . Our pre- 
liminary results indeed show so. How much loq can shift 
depends on the magnitude of V ac and co ac . The detail 
analysis in this regime will be our future direction in this 
project. 
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FIG. 5. The limit cycles in phase plane V5 — V38 under an 
extra weak ac bias Vac cos(oj ac i) with to a c being loo/3 (solid 
line) and 2.5wo (dot line). The inset on the left (right) is the 
curve of current vs. time for u) ac = u>o/3 (=2.5wo)- All other 
parameters are the same as those for Fig. 4. The current 
frequency equals to u) ac (uj ac /5) for uj ac — loq/3 (=2.5o>o)- 

In summary, we show that the trajectories of SLs in 
their phase spaces are closed curves when TCSOs occur. 
These closed curves do not depend on the initial condi- 
tions. In other words, the closed trajectories are isolated. 
Thus we show numerically the existence of the limit cy- 
cles, and TCSOs can be understood as the manifestations 
of limit cycles. Just like many other non-linear dynami- 
cal systems, TCSOs are governed by the properties of the 
unstable SSSs. According to this theory, the generation 
and motion of an EFD boundary are also the properties 
of unstable SSSs. An EFD boundary does not necessar- 
ily need to start from the emitter. It can start from an 



interior well of a SL, and it then oscillates inside the SL. 
They are universal in the sense that they do not depend 
on the initial conditions. This universal property may 
not be so obvious in previous theories [5]. We have also 
investigated the effects of a small extra external ac bias 
on TCSOs. We find that the tunneling current will oscil- 
late periodically when the ac bias frequency Lo ac is com- 
mensurate with the system intrinsic frequency loq. The 
current frequency equals either u ac or Lj ac /n, where n is 
an integer, showing periodic doubling which is a general 
phenomenon in non linear dynamical systems. 
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